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Abstract 

^: 

O , After a brief review of the BFKL approach to Regge processes in QCD and in 

supersymmetric (SUSY) gauge theories we propose a strategy for calculating the next- 
to-next-to-leading order corrections to the BFKL kernel. They can be obtained in 
terms of various cross-sections for Reggeized gluon interactions. The corresponding 
I amplitudes can be calculated in the framework of the effective action for high energy 

■ scattering. In the case oi N = A SUSY it is also possible to use the Bern-Dixon-Smirnov 

. (BDS) ansatz. For this purpose the analytic properties of the BDS amplitudes at high 

I energies are investigated, in order to verify their self-consistency. It is found that, for 

' the number of external particles being larger than five, these amplitudes, beyond one 

, loop, are not in agreement with the BFKL approach which predicts the existence of 

. . I Regge cuts in some physical channels. 

> ■ 

H 1 Introduction 

The elastic scattering amplitude in QCD at high energies for particles with color indices 
A,B and helicities Xa,^b in the leading logarithmic approximation (LLA) has the 
Regge form 

A2^2 = 2g6x^x^,nA —9^3' h^x^, , t = -<f. (1) 

The gluon Regge trajectory, j{t) = 1 +u){t), reads 

-^(-f) = i^-^r I d'~'^k « - a f In 4 - -1 , (2) 



(27r)2 ^ J k'^{q-kf 

where we have introduced dimensional regularization with D = A — 2e and the renor- 
malization point ^ for the 't Hooft coupling constant 

a = ^ (4vre--)^ . (3) 
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Figure 1: Multi-Regge kinematics 

The gluon trajectory is also known in the next-to-leading approximation in QCD [2j 
and in SUSY gauge models [3] . 

In LLA gluons are produced in the multi-Regge kinematics (see Fig. [1]). In this 
kinematics the gluon production amplitude in LLA has the factorized form 



where 

s = {pA+ PBf Sr = {kr + K-if > ql , kj. = qr+1 - qr . (5) 

The matrices are the generators of the SU (Nc) gauge group in the adjoint repre- 
sentation and C^{q2,qi) are the effective Reggeon-Reggeon-gluon vertices. In the case 
when the polarization vector e^{ki) corresponds to a produced gluon with a definite 
helicity one can obtain [4j 

C^C,iq2,qi)e;{k^) = V2^-^, (6) 

where the complex notation q = qx + iq-y for the two-dimensional transverse vectors 
has been used. 

The elastic scattering amplitude with the vacuum quantum numbers in the t- 
channel can be calculated with the use of s-channel unitarity In this approach 
the Pomeron appears as a composite state of two Reggeized gluons. It is also conve- 
nient to use transverse coordinates in a complex form together with their canonically 
conjugated momenta as 

* • ■ 9 ^ . d 

Pk = xk + lyk , Pk = xk- lyk , Pk = iTT- ; Pk = ^TTT ■ [') 

dpk 

In the coordinate representation the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 
for the Pomeron wave function can then be written as follows [1] 

^^-(^1,^2) = ^i2*(pi,P2) , A = -^^min^, (8) 

ZTT 
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where A is the Pomeron intercept. The BFKL Hamiltonian has the rather simple 
representation [5] 

Hi2 = In |piP2p + ^-rO" |Pi2p)PiP2 + (In \Pi2\^)PiP2 - 4V'(1) , (9) 
P1P2 P1P2 

with P12 = pi — P2- The kinetic energy is proportional to the sum of two gluon Regge 
trajectories u}{—\p\f) (i = 1,2). The potential energy ~ In |/?i2p is related to the 
product of two gluon production vertices C^. This Hamiltonian is invariant under the 
Mobius transformation [6] 

apk + b 

Pk — ^ , (10) 

where a, 6, c and d are complex numbers. The eigenvalues of the corresponding Casimir 
operators are expressed in terms of the conformal weights 

1 n _ 1 n 

m = + m = -+iv - - (11) 

for the unitary principal series representation of SL{2,C), with being real and n 
integer. 

The BFKL Hamiltonian can be iterated in the s-channel to account for the exchange 
of an arbitrary number of Reggeized gluons. This iteration is described by the Bartels- 
Kwiecinski-Praszalowicz (BKP) equation [7] for the n-gluon colorless composite state. 
In the Nc ^ 00 limit the Hamiltonian has the property of holomorphic separability f8] 
in the form 

H = ^J2 ^k,k+i = ^{h + h*), [h, h*]=0. (12) 

k 

The holomorphic Hamiltonian can be written as 

^ = Yl hk,k+i , hi2 = ln(piP2) + — (In P12) Pi + — (In Pi2)p2 - 2^(1) , (13) 

where iIj{x) = (Inr(x))'. Consequently, the wave function ^ fulfills holomorphic fac- 
torization [8j and there exists the remarkable duality symmetry under the transforma- 
tion [9] 

Pi pi^i+i pi+i . (14) 

Moreover, in the holomorphic and anti-holomorphic sectors, there are integrals of mo- 
tion commuting among themselves and with h jS] [lO] : 

Qr= Pklk2Pk2k3--PkrklPklPk2--PkrA'lr,h]=0. (15) 

kl<k2<...<kr 

The integrability of BFKL dynamics was demonstrated in [10] and it is related to the 
fact that h coincides with the local Hamiltonian of the Heisenberg spin model [11]. 
In the LLA the Pomeron intercept is A = 4 —Nc ln2 > and the Froissart bound 
at < cln^ s for the total cross-section ~ s'^ is violated [I]. 

In the next-to-leading logarithmic approximation the integral kernel for the BFKL 
equation was constructed in Refs. [3, 12j. Due to its Mobius invariance, the solution of 
the BFKL equation can be classified by the anomalous dimension 7 = ^ + iz^ of twist-2 
operators and the conformal spin \n\, which coincides with the number of transverse 
indices of the local operators O^^,,,^^. 

The eigenvalue of the BFKL kernel in the next-to-leading approximation has the 
form [12] 



\ti\\ ( \ti\ 
2V'(1) - V'lT+Yj-^U-^+V 



+ 4a2A(n,7), « = • (16) 
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In QCD, the next-to-leading contribution A(n,7) is a non-analytic function of the 
conformal spin |n| because it contains some terms depending on the Kronecker symbols 
5nfl and 6n,2- However, in = 4 SUSY this dependence is cancelled and we obtain 
the following hermitially separable expression [U [13] 

A(„.,) = m) + m-) - ^HilM^n , M = 7 + M , (17) 

Za/u 2 



p{M) = P'iM) + lc(2) , /3'iz) = \ 



(18) 



where 

4>{M) = 3C(3) + / (M) - 2$(M) + 2/3' (M) (^(1) - V(M)) , (19) 

and 

^ ' k + M f^^k + M y ^ ' k + M y ^ ' 

Very importantly, all these contributions have the property of maximal transcendental- 
ity [13]. The behaviour of the 4-point Green function corresponding to this NLO kernel 
in iV = 4 SUSY was investigated in [H]. The NLO conformal spins affect azimuthal 
angle decorrelations in jet physics as it was originally suggested in [15j . 

In a different context, the one- loop anomalous dimension matrix for twist-2 op- 
erators in = 4 SUSY can be easily calculated since it is completely fixed by su- 
perconformal invariance. Its eigenvalue is proportional to tp{l) — — 1), which is 
related to the integrability of the evolution equation for the quasi-partonic operators 
in this model [16j . The integrability of = 4 SUSY has also been established for other 
operators and in higher loops [iTl [18] . 

The maximal transcendentality principle suggested in Ref. [13] made it possible to 
extract the universal anomalous dimension up to three loops in A'^ = 4 SUSY [19[ [20] 
from the QCD results [21]. This principle was also helpful for finding closed inte- 
gral equations for the cusp anomalous dimension in this model [221 123] based on the 
AdS/CFT correspondence [24[ [25l 126] . In the framework of the asymptotic Bethe 
ansatz the maximal transcendentality principle helped to fix the anomalous dimension 
at four loops [27]. However, the obtained results contradict the predictions stemming 
from the BFKL equation [3l [13] . The origin of this discrepancy is related to the onset 
of wrapping effects [27]. In this framework it is, therefore, crucial to obtain more infor- 
mation from the BFKL side through the calculation of its higher order corrections to 
the integral kernel. We would like to point out that the intercept of the BFKL Pomeron 
at large 't Hooft coupling constant in A^ = 4 SUSY was found in Refs. [20[ [28]. 

In the present paper we want to formulate a program to calculate the three loop 
corrections to the BFKL kernel in the 't Hooft coupling. Our approach is based on 
the use of the high energy effective action developed in [291 [2Q] for the construction 
of the various Reggeized gluon couplings, and on the BDS ansatz [31] for scattering 
amplitudes in the A'^ = 4 super Yang-Mills theory. We begin with a short review of 
the effective action (section 2) and then turn to an analysis of the BDS formula in the 
Regge limit for the amplitudes up to six external gluons (section 3). An interpretation 
based upon known results of the high energy limit of scattering amplitudes in QCD is 
given in section 4. An outlook is presented in the concluding section. Some details of 
the calculations are presented in several appendices. 
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2 Effective action for Reggeized gluons 

Initially calculations of scattering amplitudes in Regge kinematics were performed by 
an iterative method based on analyticity, unitarity and renormalizability of the the- 
ory m. The s-channel unitarity was incorporated partly in the form of bootstrap 
equations for the amplitudes generated by Reggeized gluons exchange. But later it 
turned out that for this purpose one can also use an effective field theory for Reggeized 
gluons [29| [30]. 

We shall write below the effective action valid at high energies for interactions of 
particles inside each cluster having their rapidities y in a certain interval 

= ;^ In 1^1 , |j/ - 7/o| < ?7 , T/ << In g ■ (21) 
2 ek- \k\ 

The corresponding gluon and quark fields are 

v^{x) = -iT'^vlix) , i^ix) , ^(x) , [T^ T'] = ifabcT" . (22) 

In the case of the supersymmetric models one can take into account also the fermion 
and scalar fields with known Yang-Mills and Yukawa interactions. Let us introduce 
now the fields describing the production and annihilation of Reggeized gluons |29] : 

A±{x) = -iT''Al{x). (23) 

Under the global color group rotations the fields are transformed in the standard way 

6Vf,{x) = [Vf,{x),x], S^pix) = -xi^ix), 6A{x) = [A{x),x], (24) 

but under the local gauge transformations with x(^) ^ at x ^ oo we have 

5v^{x) = -[D^,x{x)], 6iP{x) = -x{x)iPix), M±(x) = 0. (25) 

g 

In quasi-multi-Regge kinematics particles are produced in groups (clusters) with 
fixed masses. These groups have significantly different rapidities corresponding to the 
multi-Regge asymptotics. In this case one obtains the following kinematical constraint 
on the reggeon fields 

d^A±{x) = 0, d± = n'^d^, (26) 

n'j. = (5q ± . For QCD the corresponding effective action local in the rapidity y has 
the form [29] 

S = j (fx (Lo + Und) , (27) 
where Lq is the usual Yang-Mills Lagrangian 

and the induced contribution is given by 

L,„, = Tr (Ll, + LfJ!) , Lt, = -d^A^d^A-^ . (29) 
Here the Reggeon-gluon interaction can be presented in terms of Wilson P-exponents 



Lo = itPDiP + -Tr Gl„ D^ = d^ + gv^, G^, = -{Df,, D,] (28) 



Lfn^ = -\Pewi^-9l£^v+{x')d{xY^ dlA^ 

Pexp (-g^ £^ v^ix')dix')+) d^A^ 



1 2 1 1 \ o2 

+ [v-- gv^-^v^ + g\-^v^^v- - ... ] d'^A+ . (30) 
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One can formulate the Feynman rules directly in momentum space [30]. For this 
purpose it is needed to take into account the gluon momentum conservation for induced 
vertices 

ko + kf + ... + k^ = 0. (31) 
Some simple examples of induced Reggeon-gluon vertices are 

Kit = 'f^Sa.cin+r , KTa\t = f±Kao {n^Y'^{n+r\ (32) 



rna rpc rpa rpc 

\vovii'2+ _ ^2 ( n + ^'^O ( -n + \^^ ( n + ^^i { ^20(1 aia I 02 01 apa \ (o'y\ 
^aoaia2C -Q±{n ) [U ) [Tl ) | ,+,+ ^ + I ' K'i'i) 



In the general case these vertices factorize in the form 

K::^::Zt = i-^rf± flin+r 2Tr {T^Ga,a,...a^ , (34) 

where T'^ are the color generators in the fundamental representation. In more detail, 
Gaoai...ar Can be written as [5U] 

Gaoa\...ar — , + /, + , +n /, + _i_ /,+ i , j,+ \' ' (3^) 

{io,ii,.-,v} ^'^^^O + '^h>-y'^io + '^ii + ••• + K-i> 

These vertices satisfy the following recurrent relations (Ward identities) [29j 



1.+ /\VQf\---"r+ (U+ \ 

"-r '-^aoax...arC I'^O ' "•' "-r ) 
r-1 

= —{n^Y'' ^ ifaaraiKo^Ya^illaa^+i...ar-ic{kQ-, kf_i, kf + k^ , ^itfi, A:,._i).(36) 
1=0 

With the use of this effective theory one can calculate the tree amplitude for the 
production of a cluster of three gluons, or a gluon and a pair of fermions or scalar 
particles (in the case of an extended supersymmetric model) in the collision of two 
Reggeized gluons [30] (see also earlier calculations of this amplitude in [32]). The 
square of the amplitude for three particle production integrated over the momenta of 
these particles is the main ingredient to construct the corresponding contribution to the 
BFKL kernel in the next-to-next-to-leading approximation using the methods of |33] . 
One can go to the helicity basis of produced gluons or fermions |34j . In principle it is 
also possible to calculate the loop corrections to the above Reggeon-particle vertices 
with the use of the effective action, however, in the present paper, we will use for this 
purpose the results for = 4 SUSY amplitudes presented by Bern, Dixon and Smirnov 
in [ST]. 



3 BDS amplitudes in multi-Regge kinematics 

As we have already remarked in the previous section, to find the next-to-next-to- 
leading corrections to the BFKL kernel in = 4 SUSY we need to calculate, apart 
from the amplitude for the transition of two Reggeized gluons to three particles, also 
the three loop correction to the gluon Regge trajectory, the two loop correction to the 
Reggeon-Reggeon-gluon vertex, and the one loop correction to the amplitude for the 
transition of two Reggeized gluons to two gluons or their superpartners. In this section 
we consider, as a first step, the corrections to the Regge trajectory and corrections to 
the Reggeon-Reggeon-gluon vertex (valid up to one loop) which can be obtained from 
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the multi-Regge asymptotics of the amphtude with the maximal hehcity violation, 
calculated by Bern, Dixon and Smirnov (BDS) [31]. We also investigate the six point 
amplitudes 2 — > 4 and 3 — > 3 in multi-Regge kinematics, thus preparing the comparison 
with QCD calculations to be carried out in the following section. 

The BDS formula determines the logarithm of the scattering amplitude (to be 
more precise: after the Born amplitude has been removed). In our analysis of the BDS 
formula we will, thoughout our paper, restrict ourselves, in the logarithm of the ampli- 
tudes, to those terms which are singular or constant in e, i.e. we do not (yet) consider 
corrections of order e or ^ in the logarithm of the amplitude. As a consequence, all 
results for the scattering amplitude are correct up to relative corrections of the order 
e, i.e all results should be multiplied by a factor of the form (1 -|- 0(e)). Details of our 
analysis of the BDS formula are outlined in several appendices. 




Figure 2: Elastic scattering in the Regge asymptotics 



According to Ref. ^31j, in the case of maximal helicity violation the amplitude A.^ 
with n legs in the large- A^c limit is factorized in the product of the tree result (including 
the corresponding color structure) and the simple scalar quantity M„. In the Regge 
limit s ^ (— t) the expression for M4 having the singularities in s and t-channels can 
be simplified as follows (see Figl2] and Appendix A)Q| 

M2^2 = m (^-^j r(t) • (l + 0(6)), (37) 



where fi^ is the renormalization point. 



-(t) = -^ln4+r^f^+/?(a)), (38) 



4 /x2 Jo a' V 4e 

is the all-order gluon Regge trajectory, as obtained from the BDS formula [MIEZ] (for 
a verification by comparison with explicit calculations see discussions below), and 

fj,"^ Jo a' \ 8e 2 / 2 



''da' a f-/(a') p(a') , 
— 1"- (^ + ^+'^(«))' (39) 



a a 



^As we have said before, the factor (1 + 0(6)) on the rhs is present in all our results for scattering 
amplitudes, and it will be omitted in the following. For example, for the calculation of the vertex function 
r, from the BDS formula, our neglect of order-e corrections in the logarithm has consequences: in order to 
determine the vertex function beyond the one loop approximation, it is necessay to compute, in the logarithm, 
also the higher order terms. However, such a computation is not the aim of this paper, and we will not 
go beyond the one loop approximations for the logarithm of vertex functions, Inr(t) or hiT{t2,ti,liiK,) (see 
below) . We will have to come back to this question when calculating the full higher order corrections to the 
BFKL equation. We thank V. Del Duca for discussions on this point. 
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is the vertex for the couphng of the Reggeized gluon to the external particles. 
The 't Hooft coupling is defined as in eq.(3): 



27r 



and the small parameter e is related to the dimensional regularization 4 
cusp anomalous dimension 7(a) is known to all loops |201 [35} [23] 

7(a) = 4a - 4(2 a^ + 22C4 a^ + ... , 
and the functions /3(a), 5{a) and C(a) read [31] 

(5{a) = -C3a2 + (6C5 + 5C2C3)a' + -, 



(40) 

4-2e. The 

(41) 



5{a) 
C{a) 



- C4 + . 
2 



where C('^) is the Riemann ^-function 



C(n) 



A:=l 



(42) 



(43) 



Written as in Eq. (j37p we can see that the asymptotic behavior of the M2^2 BDS 
amplitude corresponds to the Regge ansatz with the gluon trajectory j = 1 + uj{t) 
given by the perturbative expansion 



1 1 
In ^ + - 



+ 
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C4 In 



a + 



C2 



In —TT 

2 2e 



1 

3e 



+ 



fi^ ze / 2 
6C5 + 5C2C3 



a-^ + .... 



(44) 



The first two terms in this expansion are in agreement with the predictions in Refs. [U 
[3]. Note that in Ref. [3], where the BFKL kernel at NLO was calculated in = 4 
SUSY, initially the MS-scheme was used, and only later, in Ref. [13] the final result was 
also presented in the dimensional reduction scheme (DRED). The NLO terms in the 
BDS expression for uj{t) can be obtained from Ref. [3j by converting it to the DRED 
scheme, where, apart from the finite renormalization of the coupling constant, one 
should also take into account in the loop the additional number 2e of scalar particles (for 
details see Appendix A and the recent paper [38j jl. The 0{a^) term in a;(t), extracted 
from the BDS amplitude [36]), corresponds to the three-loop correction to the gluon 
Regge trajectory needed when calculating the next-to-next-to- leading corrections to 
the BFKL kernel in this model. Strictly speaking the Regge asymptotics of scattering 
amplitudes corresponds to a different order of taking two limits e — > and s ^ 00, but 
it is probable that they can be interchanged. 

It is noteworthy to point out that the expression for M2^2, derived in the Regge 
kinematics, in fact, is valid also outside the Regge limit. That is to say that, when 
analysing the BDS formula for the logarithm of the 2^2 amplitude (see Appendix 
A), we do not make use of the high energy limit. In particular, the amplitude can also 
be written in the dual form 



2-+2 



r(.). 



(45) 



^ We thank A. V. Kotikov and E. M. Levin for helpful discussions regarding these redefinitions. 
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Figure 3: Production amplitude in the multi-Regge regime 



After having described the elastic scattering amplitude we now focus on the BDS 
production amplitude (see Fig. [3|). The analysis of lnM2^3, described in Appendix 
B shows that the amplitude (as before, up to the corrections (1 + 0{e)): see also the 
footnote on p. 7) can be written in the factorized form 

M2^3 = r(ti) (^-^j T{t2,h,ln-K) i^-^j T{t2), (46) 

with 

i-s) 

Since Eq. ()46p is exact it is also valid in the multi-Regge kinematics 

- s > -si, -S2 > -ti ~ -t2 ~ -K , (48) 

where all invariants s,si,S2,ti and t2 are negative. Due to the correct factorization 
properties of this amplitude the Reggeon-particle-particle vertex T{t) in Eq. (f46l) is 
exactly the same as in the elastic amplitude in Eq. ([39|) . The gluon Regge trajectory 
uj{t) in Eq. (I46p also coincides with the one discussed above. The new component is 
the Reggeon-Reggeon-particle vertex. Its logarithm is given by 

1 T.. IN 7(a) , 2 1 r da' a fj(a') (5{a') ^, 

Inr t2,ti,ln-K = -^In^ — -- In - + + 5 a' 

lb //^ 2 Jo a' a' \ 4e^ 



da' 






In 4 




a' 


F 


da' / 


Jo 


V [ 



- 4? '"^ ^ - ^c. - i (.(.)..<.)-/; ^ (^ . m) ) .n ^ . m 

It is now possible to analytically continue this 2 — > 3 production amplitude to the 
physical region where the invariants s, si and S2 are positive (see Fig. H^) 

M2^3 = r(ti) (^^^ j r(t2,ti,inK-i7r) (^^l^j r(t2). (so) 

A similar continuation to another physical region can be performed in the case 
when s is positive but si and S2 are negative (see Fig. Hja): 

M2^3 = r(ti) (^-^j r(t2,ti,in/t + i7r) (^-/j r(t2). (si) 
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Pa' Pa 



k k 



Pb' Pb 




Figure 4: Physical channels for the one particle production amplitude 






+ C2 



Figure 5: Dispersion representation for M2^3, exhibiting its analytic structure 



Using a 'dispersive' representation illustrated in Fig. [5] (in the following, we will 
refer to this type of representation as 'analytic' representation), it can be easily verified 
that in all physical regions (including the crossing channels with s,si < 0,S2 > and 
s, S2 < 0, si > 0) the amplitude can be written as follows 



M2- 



-si 



-SK 



W(i2) 



with the real coefficients Cj 

ci = \T{t2,ti,ln-K)\ 

C2 = |r(t2,ii,ln-K;)| 

where cpr is the phase of F, i.e. 



ci + 



-S2 



sin7r(w(ti) — (pr) 
sin7r(a;(ti) — ^(^2)) 

sin7r(u;(t2) - 4>r) 
sin7r(w(t2) — ^{ti)) 



F(t2, ii, In K — in) = |F(t2, ^i, In —k)\ e 



C2(52) 

(53) 
(54) 

(55) 



In this dispersion-type representation for all physical channels we use the on-mass shell 
constraint for the produced gluon momentum 



K = k'i = {qi- q2f , 



(56) 



where k±_ is its transverse component {k±^pA = k±pB = 0). In this case the amplitude 
M2-3 does not have simultaneous singularities in the overlapping channels {si,S2), in 
an agreement with the condition of the gluon stability (this will be discussed further 
in section 4.2). 
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The fact that there exists a solution for the coefficients ci and C2 proves that the 
scattering amphtude, derived from the BDS formula, has the correct analytic structure 
in all physical regions. In particular, it satisfies the Steinmann relations (a somewhat 
more detailed discussion of analyticity properties will be presented in the following sec- 
tion). We therefore conclude that the BDS amplitude for the production of one particle 
in multi-Regge kinematics has the correct multi-Regge form. This is encouraging, and 
we proceed now to study the production of two particles in multi-Regge kinematics, 
for which we use the M2^4 BDS scattering amplitudes. 




Figure 6: Production of two particles 

We have first checked that the planar BDS amplitude for two particle production 
having singularities only at positive values of the invariants 5,51,52,53,^1,^2,^3 has 
the correct multi-Regge form in the multi-Regge kinematics in the region where all 
invariants 5,51,52,53,^1,^2,^3 are negative (see Fig. [6] and Appendix C) 

r{t2,ti,ln-Ki2) { ] r(t3,t2,ln-K23) ^- (57) 



and the quantities 

_ (-5l)(-52) _ (-52)(-53) , . 

— K12 — , — K23 — • lOOj 

— 5012 —5123 

are fixed together with ti. The invariants S012 and 5123 are the squared masses of the 
corresponding three particles in their center-of-mass system. 

In a similar way to what we did in the 2 — > 3 case, the BDS 2^4 amplitude can 
be analytically continued to several physical channels, each of them corresponding to 
different signs of the invariants s, soi2, 5123, 5i, S2, 53. To begin with, in the region (see 

5, S012, 5123, Si, 52, 53 > (59) 

the amplitude has the form 

M2^4 

r(ti)r(t3) " 

-j^j r(t2,ti,lnKi2 - ivr) T{t3,t2,lnK23- itt) y-^j ,(60) 

where we can replace K12 and K23 by their values on the mass shell 

5152 -2 S253 ,-2 

Ki2 = fci^ , K23 = fc2± • (61) 

5012 S123 
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Figure 7: Physical channels for two particle production 



Here ki = qi — q2 and k2 = q2 — Qi are the momenta of the produced particles (see 
Fig.Ej). 

In the physical region, represented in Fig. [Tja, where 



S, Sou, S3 > , Si, S2, Sl23 < 0, 



(62) 



one can obtain, for the multi-Regge asymptotics of the BDS amplitude, the following 
expression 



M2- 



-si 



T(t2,ti,lnKi2 + ivr) , „ 



S2 



T{t3, t2,lnK23 - ivr) 



-S3 



In a similar way, in the region (see Fig. [7]:) 

S, Sl23, Sl > , S3, S2, Soi2 < 

we have 

M2^4 



(63) 



(64) 



-si 



r(t2,ti,ln«;i2 - ivr) 



-S2 



T{t3,t2,ln K23 + i-n-) 



-S3 



.(65) 



We can now attempt to write M2^4 in all physical regions in terms of the following 
dispersion ansatz (represented in Fig. [8|) 

M2^4 ^-S012^^,"(*2)-"(*^) ( -SKi2K23 y^''^ 



r(ti)r(t3) 
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Figure 8: Analytic representation of the amplitude M2^4 



— S2 \ ^ y -^f / —5123^23 \ ^ 'I —5/^12^23 \ 



d3 

(66) 



with the real coefficients dj=i,2,3,4,5- Here 

Ki2 = (^L - qhf , «^23 = (9*2 - qsf ■ (67) 

By comparing the above 'dispersive' representation with the previous expressions 
for the BDS amplitude in the three physical regions Figs[7^,b,c it is possible to extract 
the coefficients. They read 

di = Ci(ki2) Ci(k23) , 
d2 = C2(ki2) C2(k23) , 
d3 + d4 = C2(ki2) Ci(k23) , 
d5 = Ci(ki2) C2(k23) 

where, in fact, ci{k) and C2(k) were defined in Eqs. ([53]) . ([5 



^ / N IP I sin7r(a;(ti) - ^na) 

Ci[Ki2) = |ri2|^ nT~\ TTTT' (69) 



sm7r(a;(t2) — i^(i3)) 

C2(K23) = |r23| """/"^?"'^H ^ (72) 



sin7r(a;(ti) 


-u;{t2)) 


sin7r(a;(t2) 


- (PT12) 


sin7r(a;(t2) 


-w(tl)) 


sin7r(Lij(t2) 


- '/Tag) 


sin7r(a;(t2) 


-u;(t3)) 


sin7r(ti;(t3) 


- <Pr23) 


sin7r((x;(t3) 


-u;it2)) 
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with (cf.dMD) 

= r(t2,il,ln-Ki2) , = r(t3,t2,ln-K23) . (73) 

We note that for the coefficients da and d4 only their sum can be determined from 
the three physical regions previously discussed. However, an attempt to fix separately 
these two coefficients from the multi-Regge asymptotics in the physical region (see Fig. 

m) 

S, S2 > , Si, S3, Soi2, Sl23 < 

leads to a disaster: the corresponding equations do not have any solution. The reason 
for this is that the BDS amplitude in this region does not have the correct Regge 
factorization (see the discussion in section 4). According to Appendix C its asymptotics 
here is 



M2- 



r(ti)r(t3 

'-si 



C 



11^ 



iir] 



T{t2,ti,ln Ki2 

where the coefficient C is given by 



-S2 



r(t3,t2,ln K23 



ITT) 



-S3 



c 



exp 



ITT 



In 



{k 



1 + A;2)V^ 



,(74) 



(75) 



The fact that, for this region, we find no solution for the coefficients di indicates that, 
in this region, the BDS amplitude does not have the correct analytic structure. In 
section 4 we will show, by comparing with explicit calculations of the high energy limit 
of scattering amplitudes, that in the BDS formula a piece is missing. This piece belongs 
to a Regge cut singularity, which - apart from the one-loop approximation - does not 
fit into the simple exponentiation of the BDS ansatz. In Appendix C we write down 
the amplitude M2^i also in the quasi-multi-Regge kinematics, where the variable S2 
is fixed. 




Figure 9: Three particle transition 

To continue our analysis of the BDS six point amplitude we now discuss the asymp- 
totics of M3^3 (see Fig. [9]) . According to Appendix D in the multi-Regge region where 
all invariants s, si, S3, S13, S02, S2 = ^2 are large and negative, its asymptotics is similar 
to the corresponding asymptotics of the M2^4 amplitude, i. e. 

r(ti)r(t3) 

/-Si\-(*l) /_s2\-(*2) /_s3X-(t3) 

y~l]TJ T{t2,ti,ln-Ki2) y-j^ j T{t3,t2,ln-K23) \^-j^J . (76) 
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This BDS amplitude can be now analytically continued to the physical region where 





Figure 10: Physical regions for the amphtude M; 



3^3 



the invariants 5,51,53,512,502,^2 are positive (see Fig. [TOk). The resulting amplitude 
can be written as 



r(ti)r(t3) 

j r(t2,ii,lnKi2 - ivr) 2-j r(t3,t2,lnK23 - ^v^) l^-^j ,(77) 

Similarly, the analytic continuation to the region where 5i, 512, t2 < and 5, 53, 502 > 
(see Fig. llOb) is of the form 

M3_3 



r(ti)r(t3) 

-5l\'^'^*^^ /-S2\'^''*^^ /-53\'^*-*^^ 

-^j r(t2,ii,lnKi2 + ivr) r(t3,t2,lnK23 - l^-^j . (78) 

Finally, the continuation to the region where 53,502,^2 < 5,5i,5i2 > (see Fig. 
[TOh ) reads 

M3^3 



r(ti)r(t3) 

-5l\'^(*i) /_S2\'^(*2) /-53\'^^*^^ 

— 2-j T{t2,ti,lnK.i2 - in) y—^J T{t3,t2,ln K23 + iir) {^—^ j -(79) 

As it was done in the M2^4 case one can write the dispersion relation for M3^3 
valid in these physical regions, which includes five contributions, shown in Fig. [TTl 
and calculate the real coefficients 61,62,63 + 64,65. But once again, we find that it 
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Figure 11: Analytic representation of the amplitude M; 



3^3 



is impossible to fix separately the coefficients 63 and 64 from the BDS amplitude, 
calculated in the physical region where si, S3, S13, S02 < and s, ^2 > (see Fig. [TOli) 



M3- 



r(ti)r(t3) 



C 



-si 



r(t2,ti,lnKi2 + ivr) 



-S2 



aj(t2) 



r(t2,ti,ln K23 + i7r) 



-S3 
/.2 



where the phase factor C' is 
C = exp 



lK{a) {qi- (92 - qzf 
y—i'K) m ■ 



(gi + 93 - 9*2)^ 



,(80) 



(81) 



The reason for this drawback is the same as before: the absence of a correct Regge 
factorization for the BDS amplitude. In the next section, using the BFKL approach, 
we shall discuss the reason for this problem. Namely, the BDS amplitude does not 
contain the Mandelstam cut contributions plotted in Fig. [121 



4 Regge cuts and breakdown of factorization 
4.1 Regge pole models 

The results of the previous section can best be understood if we confront them with 
the known high energy behavior of QCD scattering amplitudes in the Regge limit. In 
the LLA, the high energy behavior of the QCD scattering amplitudes is the same as 
in the supersymmetric case. We proceed in three steps: we first review the findings 
for models containing only Regge poles. We then summarize the results obtained in 
gauge theories, and finally compare with the scattering amplitude derived from the 
BDS formula. 
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A key element in analyzing the high energy limit are the Steinmann relations [39] 
which forbid the existence of simultaneous energy discontinuities in overlapping chan- 
nels. As an illustrative example of the Steinmann relations, consider the 2 — > 3 am- 
plitude shown in Fig. [3 obviously, the produced particle in the central region can 
form resonance states with particle A' or with particle B', but not simultaneously with 
both of them. As a result, in the physical region the scattering amplitude cannot have 
simultaneous discontinuities in the energy variables si and S2- The way in which this 
restriction is implemented into scattering amplitudes is that, in the double Regge limit, 
the signatured amplitude can be written as a sum of two pieces, one of them with cuts 
in the si and in the s channels, the other one in the S2 and in the s channels. In general, 
there are cut singularities both in the right and left half energy planes, and one has 
to form signatured combinations. Decompositions of this kind have first been derived 
from simple models which contain only Regge poles (massive scalar tp^ theory [30] or 
the dual Veneziano 6 point amplitude, Bq [S]), and from studies of dispersion relations 
and generalized Froissart-Gribov partial wave representations [32]. For 2 ^ 4 or 3 ^ 3 
amplitudes, we have five independent terms, and for scattering processes with higher 
number of legs the number of terms grows rapidly. 

For models which contain only Regge poles the general structure of the signatured 
2 — > 3 amplitude is: 



A 



2^3 



pA{h)PB{t2) 



+ 



+ ^1^2 -2 



-£2 
/x2 



a{t2)-a{ti) 



+ nT2 



a(t2)-a(il)' 



a(t2) / s \ 



Vi{ti,t2, k) 
t^2(tl,t2,K).(82) 



Here a{ti) denotes the trajectory function of the Regge pole in the ti exchange channel, 
Ti (T2) are the signatures of the ti {t2) channels, and as usual, (— s)° = (|s|)"e~*™. 
In this representation, the energy singularities are explicit, i.e. all phase factors are 
contained in the energy factors, and the functions Vi are real valued functions. 

With the abreviations 

ai = a{ti),aij = a{ti) - a{tj) (83) 
and with the signature factors 



e + TiTj 



(84) 



we can rewrite the expression for ^2^3 
^2^3 



pA{tl)f3B{t2 



+ 



f\s\Y\ , Vlih,t2,K) 

fJ- J VA* / sm7rai2 
fi^ J \lJi'^ J sm7ra2i 



(85) 



where the vertex function Vi is proportional to Vi . The generalization to the signatured 
2^4 amplitude (consisting of five different pieces) is illustrated in Fig. [HI and from 
Eq. ([66]) one easily obtains the analogue of (l85l) |41j : 



A 



2-*4 



PA{tl)fiB{h) 

SlW'^ f\soi2\V f\s\Y\ , , Vi{h,t2,Ku)Vi{t2,h,K23) 

2 I ?12?23?3 



sm 7rQi2 sm 7rQ23 
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fj,'^ J \ fj,'^ J KfJ'^J '"^^"'^^''i sin 7ra32 sin 7ra2i 

+ (My /^kody flsiy ^2 (tl , ^2 , ) Vl {t2 , h , ^23 ) 

\ fi^ J \ jj? ) KfJ-'^J sin7ra2 sin 7ra2i sin vrais 

+ (My (Ifl^y (\s\y\^^,^^,^ Sm7TasV2{h,t2,Ki2)Vlit2,t3,K23) 

\ ^J? ) \ fx'^ J \fJ.^J sin7ra2 sin 7ra23 sin vrasi 

+ (N)" 04)'- (q)"faA.fa '^-"-.'-""'^'fe''^'''-^> . (86) 

V /U^ / V / VA* / sm 7ra32 sm 7rai2 

The analogue for the 3 — > 3 process again consists of five pieces which are shown in 

Fig. [in 

From the discussions of these Regge pole models it has also been recognized that 
the analytic decomposition into a sum of terms in ()85p and (I86p is consistent with a 
factorizing form. For the 2^3 case we can write: 



A2^3 ^ /N\"^ . ..^^N2r"^ 

pA{ti)pB{t2) V 



ClVr,r,{h,t2,K) (^^^ (87) 



Here the important point to be stressed is that the new production vertex function 
FriT2 contains phases (in contrast to the real- valued functions Vi in (f85]) ). and it has 
cut singularities in the K-plane. Similarly for the 2 — > 4 case we have 



A /loIX"^! / \ c \ \ /lol\"3 



^1 ( "TT ^rir2(tl,t2,Kl2)6 — Karg (^2 , ^3 , '«23) 6 — 



/5yl(il)/3B(t3) V /"V V /"V V /" 

with the production vertex function from (j87p . As a result, for this class of Regge-pole 
models the production amplitudes, in the multi-Regge limit, can be written either in the 
'analytic' form (sum of terms with simple analytic properties and real-valued vertex 
functions Vi) or in the 'factorized' form (with the production vertices V containing 
phases and singularities in k). 

Let us comment on the planar approximation. Planar amplitudes have right hand 
cut singularities only, and in the physical region where all energies are positive, their 
phases follow from the signatured amplitudes in (I85p or ()86p by simply dropping all 
'twisted' terms containing factors Tj. One can easily verify that, in the physical region 
where all energies are positive, these planar amplitudes can also be written in the fac- 
torized form (fHTI) and (188]) (with vertex functions ^^(ti, ^2, ^12) being slightly different 
from the signatured ones, Ftit2 (*i) *2> '^12))- When analytically continuing into the un- 
physical region, where all energy variables are negative and well separated from their 
threshold singularities, all phases inside the production vertex V{ti,t2, K12) disappear, 
the vertex function turns into a real- valued function, and the factorized form remains 
valid. However, in the physical region where s,S2 > and si, S3, S012, -5123 < the 
factorized form is not valid, and the structure of the amplitude is more complicated. 



4.2 High energy behavior in Yang Mills theories 

Let us now turn to QCD. Throughout this section we will restrict ourselves to scattering 
amplitudes with odd signature in all t-channels. Compared to the Regge pole models 
discussed in the previous subsection, the situation is slightly more complicated since 
also Regge cut pieces appear in some of the t-channels. In the LLA the real part of the 
2 — > n scattering amplitude is well known to have the factorized form of Eq. and 
it is in agreement with our previous result in (j87p and (j88p . However, when turning to 
the imaginary parts (i.e. to the energy discontinuities) of the production amplitudes, 
a new piece appears which destroys the simple factorization property. The best way 
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of understanding the appearance of this new piece is the use of s-channel unitarity in 
the physical region where all energies are positive. 

Starting from the analytic representation of the scattering amplitude A2^n, it is 
possible to determine, in QCD, the partial waves from energy discontinuities and uni- 
tarity equations |43l I44j . As the simplest example, let us consider, in the LLA, the 
2 — > 3 amplitude, consisting of the two terms illustrated in Fig. \5\ Anticipating that, 
in the 2 — > 3 process, there are only Regge pole contributions, we start from the ansatz 



T{ti)T{t2) ht2 



f\s\Y\ , Vi{ti,t2,K 



2 



0^21 / \ \\UJ-i 



S2\y f\s\Y\ , V^2(tl,t2,Ac) 



a 



/i^ / VAi^/ " " sm7ru;2i 



with (cf.dM 
and 



(91 - q2Y (90) 



r(ti) = gdx^x^, , T{t2) = gSxsX^, . (91) 
Here we used ai{ti) = 1 + a;(tj) and 

uJi = uj{ti), LOij = uj{ti) - uj{tj), (92) 

and the signature factors can be written in following form: 

C^ = e-*"^' + 1, Cij = e-^""'^ + 1. (93) 



Taking the discontinuity in si, only the first term in (j89p contributes. Making use of the 
unitarity equation and invoking, for the ladder diagrams in the ti channel, the BFKL 
bootstrap condition we find the partial wave Vi in the LLA. For a definite helicity Vi 
has the form 



Vi = gTTC{q2, qi) {-{ioi- u;2) - -(In — - -)) (94) 



1 a K 1, 

-(^i-^2)--(ln-^ 

with C{q2,qi) being the production vertex from ([6]). Similarly, the discontinuity in S2 
leads to 

V2 = gnC{q2,qi) Q(^2 - ^i) - ^0"^ " ^)) • (95) 

A comment may be in place on the term In k in Vi and V2: it indicates that, in contrast 
to massive field theories where the V^'s are analytic functions near k = 0, in massless 
theories this is no longer the case. Therefore, when computing the discontinuity in si or 
S2 of ^2^3) there is, at first sight, an uncertainty in handling the cut in k. It turns out 
that the correct prescription for computing the discontinuity in si or S2 in the physical 
region, is keeping k = (qi — 5*2)^ fixed. This can be derived either from a direct analyis 
of Feynman diagrams where the Steinmann relations dzscs^dzscsj ^2^3 = are fulfilled 
explicitly. Alternatively, one can consider the massless Yang Mills theory as the zero 
mass limit of a nonabelian Higgs model where the gauge bosons are massive: before 
the zero mass limit is taken, the vertex functions are analytic near k = and there 
is no ambiguity in computing the energy discontinuities. As a result, in the physical 
region the singularities of in k of Vi and V2 are not related to singularities in si or 
S2- We also mention that, in the next-to-leading approximation, the functions Vi and 
V2 contain an additional dependence on Iuk, which, again, does not contradict the 
Steinmann relations 1331. 
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Inserting these expressions into Eq. (j89p . using ()93p . and restricting ourselves to 
the planar approximation, we find for the real part (apart from the color factor): 

'^'-i\s,\r 9Ciq2,qi)i\s2\r (96) 



in agreement with (4). In particular, the infrared singular pieces in Vi and V2 cancel. 
As a further test, one could also compute, from the corresponding unitarity equation, 
the single discontinuity in s: here both partial waves Vi and V2 contribute, and the 
result is in agreement with (j94p and (jDSp l. 

For the 2^4 amplitude we start from an ansatz which is slightly more general 
than (j86p . In order to account for the Regge cut in the t2 channel, we introduce, in 
the t2-channel, the Sommerfeld- Watson integral / dw'2/2'Ki: 

^2^4 2s f duj2 



r(ti)r(t3) tit2t3J 27ri 

1^012! 1^1 , , ^i(*i'*2,t3,^i2,'«23;^) 

?12'?2'3?3 



/i^ / \ fJ-"^ J V/U^/ ' ' ' sm 7rwi2' sm 7ra;2'3 

I \ UJS—Ujij /I I \ UjL—UJl /I I \ UJI TTT I X J. X I \ 

^ \S\2Z\\ \^\\ f t t ^^2(tl,t2,i3,'«12,K23;'^2j 

432'?2'1?1 



/i^ / V /^^ / Vm^/ " ' sm7rti;32' sm7ra;2'i 

'S2|^"^""VN012|^"'""VNI^"^ t t W^3(il,t2,t3,«:i2,«:23,'«123;'^^; 



,3 



^ ) \ ) J ' ' sm7rci;2'i smvrwia 

S2|\"^""VNl23l\"'~"VNh"' . t t 1^4(il,i2,i3,Kl2,K23,Kl23;a;^) 



+ 1^ Til ^2'36ie: 



V A*^ / V / Va* / sm 7rL<j32' sm 7ru;i2' J 

with the partial wave functions, denoted by VFi=i, 2,3,4,5 and to be determined from sin- 
gle energy discontinuity equations. The partial wave functions W3 and W4 also depend 
upon the additional variable K123 = {ki + ^2)^- We have the five single discontinuities 
in si, S2, S3, S012, and S123 which allow to find the partial waves VFi=i^2,3,4,5- In leading 
log accuracy the results are the following: 

Wi = Viih,t2,lil2)^ Vi{t2,t3,K23), (98) 

UJ2 — UJ2 

W2 = V2{ti,t2,Ki2)^ V2{t2,t3,l^23), (99) 

UJ2 — UJ2 

sinvrwi , n 1 / x 

W3 = V2{ti,t2,Kl2)— Vi{t2,t3,K23) 

sm TTL02 '^2 ~ "^2 

- sin7r(w^-a;i)(ycut-^p), (100) 
sinvrws , , 1 , , 

^^4 = -. V2(tl,t2,Ki2 ^ V^l i2,t3,«^23 

sm TXUJ2 — W2 

- sin7r(u;2 - ^3) (Kut - Vp) , (101) 

W5 = Vl{ti,t2,Ki2)^ V2{t2,t3,K23)- (102) 

UJ2 ~ '-^2 

The three amplitudes Wi, W2 and W5 are products of the production vertices Vi in 
(f94|l and (f95]) . found in the 2 — > 3 case, whereas the amplitudes W3 and W4 contain, 



•^We emphasize that the same results are obtained if one starts from the double discontinuities in s and 
Si: using unitarity conditions and making use of generalized bootstrap conditions, one again arrives at (|94p . 
This is a crucial test of the selfconsistency of this 'unitarity-based approach'. 



20 



in addition to the products of production vertices Vi, the extra pieces, V^ut — Vp which 
will be defined in the following. The term V^ut contains Regge cuts and cannot be 
written as a simple product of vertices for the two produced gluons. It takes the form 
of BFKL-like ladder diagrams in the color octet channel, and it is illustrated in Fig. 
12] (left figure): 



Vr 



cut 



t2Nc 4 



<Pk<Pk' 



(2^)6 {k + ki)' 



;C{k,k + ki] 



X G^^'^Xk, q2 - k; k', q2 - k'-uj'^)C{k' - ^2, k') — 



93 



{k' - k2f 



(103) 



Here C denotes the effective Reggeon-Reggeon-gluon vertex given in (6), and G^s-*) 

IS 

the BFKL Green's function in the color octet channel, satisfying the integral equation 



J2G^^^\k,q - k;k' ,q - k' ■J2) 
(27r)35(2)(A;- A:' 



A;2(g - kf 



+ 



(104) 



where K denotes the BFKL kernel in the color octet channel, containing both real 
emission and the gluon trajectory. In lowest order in the coupling, and for equal 
helicities of the two produced gluons, Vcut equals: 



(0) ^ 2 C'(g2,gi)C'(g3,g2) 
2^ 



uji + uj'i + a \ \n. 



^123 



(105) 



The term Vent not only violates the simple factorization of Regge pole models, but 
also, when computed beyond the one loop approximation, will be shown to disagree 
with the BDS formula. Finally, the subtraction term Vp removes the Regge pole piece 
inside Vcut, and it is of the form: 



yp = 9 



,2<^( g2,gl)C(g3,g2) 

4a;2 



wi + u;2 + o In 



1 



W2 + ^3 + o In 



K23 



^2 — ^2 
1^ 



(106) 



Before we compare with the BDS formula, let us remark on a few features of these 
leading order QCD results (for details, see Appendix E). From now on, we specialize 
on the planar approximation, i.e. in the signature factors in eq.(|93p we only retain the 
phases. Inserting the results of (|98p - ()102p into the full amplitude ()97p we can derive 
the results for different kinematic regions. 

Beginning with the physical region where all energies are positive, one finds that 
the sum of the Regge pole terms can be written in the simple factorizing form ()88p . In 
particular, the Regge cut pieces contained in and W4 cancel completely, and the 
real part of the scattering amplitude coincides with (4). 

Next, in the unphysical region where all energies are negative and all phases disap- 
pear, again, the Regge pole contributions can be written in a simple factorizing form, 
and the cut pieces in W3 and W4 cancel. 

Finally, we go into the physical region where s,S2 > and si, S3, 5012, -5123 < 0. 
Nonzero phases appear only in s and in S2- After some algebra we obtain: 



A 



2^4 



2s 



r(ti)r(t3) tit2h 



9^C{q2,q^)C{q^,q2){\s^\r{\s^\) 



1^3 
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Figure 12: BFKL contributions to the amplitudes M2^4 and M; 



3^3 



e-^'^"2(|s2ir[l + i| (^u;i+c^2 + a(^ln^- + + ^2 + a (in ^ - ^ 

-2zvr^|^(e-|.2ir^^cn..(107) 

In the last term, Vcutt it is possible to factor out the gluon trajectory (details are 
presented in [i7]): 

(108) 

where in the one loop approximation ()105p 



5^(^(92,^1)^(^3,92) / Tr^(e~'''|s2|)''2T4«t,redW 



9 



2C'(92,9l)C(g3,92) 



A '^123/^^ , 1 1 I Al^ 2, ^ 

a In — + - + 0(0^ In S2) 



(109) 



and the two loop and higher order terms of Vcut,reduced are infrared finite [17] . 

Inserting (llOSp into ()107p we see that all terms on the rhs of (|107p are proportional 
to the common phase factor e"*'^'^^; 



2^4 



2s 



(kiir(k3|)"^(N2|)"^5'C(g2,gi)C(g3,92)e 



r(ti)r(t3) tit2t3 



vr / / K12 

"2 V ^ +^2 + a yln — 



1\ I K23 

+ u;3 + 0^2 + a In ■ 



-2iTT / — ^(e 



^2\) ' cut, reduced 



(110) 



and the coefficient in the square brackets is infrared finite. This shows that, in LLA, 
the imaginary part of ^2^4 is infrared singular, but the singularities are assembled in 
the phase factor e~™^'^ . This observation will be important when comparing with the 
BDS formula. 
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A completely analogous discussion applies to the case 3 — > 3 (Figs. [U [TTj) in the 
multi-Regge region (a more detailed discussion is given in Appendix E). Again, the 
scattering amplitude consists of five terms, and two of them contain the Regge cut 
piece: 



/ 



7— (e \S2\) ^Ucut 
ITTI 



with 



X G^^^\Kq2-k-k\q2-k'-J^)C{k'-k2,k')j^^^^. (Ill) 

In lowest order (and for equal helicities of the produced gluns) this Regge cut piece 
equals: 

g'^C{q2,qi)C{q3,q2) . K12K23 
Ucut = ^77 aln —2-2. (112) 

2uj2 [qi + 93 - 12) q2 

Proceeding in the same way as for the 2 — > 4 amplitude, on derives results for the 
scattering amplitude in the different kinematic regions. In the region where all energies 
are positive, we find the same factorization as for simple Regge pole models, i.e. the 
Regge cut pieces cancel. In the region s, S2 > 0, si, S3, S13, so2 < the Regge cut piece 
appears. First we rewrite it in the same form as in (jlOSp : 

^(e-^ns2|)'^^t/cni = 5'C(<72,gi)C(g3,g2)(e-*ns2|)'^^/ ^(e-^1s2|)"^C/c.t,.edW 

(113) 

with the infrared finite one loop approximation 



9^C{q2,qi)C{q3,q2) / ^(e '^\s2\Y''^U cut, reduced 



^2 C{q2,qi)C{q3,q2) 



^^^^ ^ ^12^23 )^Oia^lns2) 

(gi +<?3 - 92)^^2/ 



(114) 



As in the case of the 2 — > 4 amplitude, the higher order corrections (denoted by 0(0^)) 
are infrared finite. With this result, the 3^3 amplitude can be written in the form 
rcf.ffTTOll): 

r(t')T{t3) = ^J^(l^iir(k2|)"^(k3|)"«5'C(g2,gi)C(g3,c?2)- 



.IT 

I + I— [loi + a 



1 '^12 l\ , , A K23 1 

In + 1^3 + a In 



-2i7r 



duj'2 



{e-'^\s2\r^Ucut 



reduced 



(115) 



Note that, in contrast to the 2 — > 4 case, there is no common infrared singular phase 
factor, e~*'^'^2, and the square bracket term on the rhs of (jll5p is infrared finite. This 
shows that the infrared structure of ^3^3 is quite different from A2_>4. 
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4.3 Comparison with the BDS formula 

Let us now return to the BDS amplitude discussed in the section 3, take the leading 
logarithmic approximation and compare with the results discussed in the previous 
subsection. In the leading logarithmic approximation we retain, in the exponent In M„, 
only the lowest order (in powers of a) of the coefficients of the energy logarithms, and 
the lowest order of the real and imaginary parts of logarithms of the vertex functions. 
In r. In the case of the 2 — > 2 scattering process the coefficient of Ins is given by 
the gluon trajectory function (eq. (f38l) ). and the leading order coefficient is the term 
proportional to a. The logarithm of the vertex function is given in (j39p : the lowest 
order term is of order a, and since t is negative, there is no imaginary part. Therefore, 
in the leading logarithmic approximation we put InF equal to zero (note that M4 
multiplies the Born approximation which contains a reggeon-particle-particle vertex of 
the order g). 

Turning to the case 2 — > 3 in the physical region, we use (I46p and (I49p (see also 
Appendix B). The new element, the logarithm of the production vertex, starts with 
terms of the order a, and the real part can be neglected (i.e. the absolute value of 
T{t2,ti,K) can be put equal to unity). But, depending upon the kinematic region, 
terms with ln(— k) may lead to imaginary parts of order a which have to be kept. In 
the region where all energies are positive the relevant terms of order a are (see ([55]): 

= ^(^(ti) + u;{t2)) + ^{}^^-\)- (116) 

In order to compare with the QCD results we use ([52]) . ([53|) . and In ([52]) we 

approximate the factors — > 1 etc, and for the real coefficients ci and C2 we 
obtain: 

l(u;(ti)-^(t2)-a(ln^-i)) 1 (^(t^) _ ^(t,) _ a(ln ^ - i)) 
= 77~\ 77~\ ' ^2 = 7—: 7—: (117) 

LO{ti) - LL>{t2) Uj{t2j - Uj{tl) 

which agrees with the leading log result in ([89]) , ([M]) , and ([95]) . In the unphysical region 
where all energies are negative we have no imaginay parts and again find complete 
agreement with the results of the previous subsection. 

In the case of 2 —> 4 we begin with the physical region where all energies are 
positive. Using eqs. ([66]) - ([68]) and applying the same arguments as for the 2 — > 3 case, 
we find 



ru ^ru \ - (e"^"^r (ci('^i2) + C2(Ki2)) (e , , 

r(ti)r(t3) [i^ fi^ 

(cl(K23) + C2(A^23))(e-^'^^)'^^ (118) 

quite in agreement with the LLA of the Regge pole part in ([E.3p . Since, in the QCD 
calculation for this kinematic region, the Regge cut pieces cancel, there is no conflict 
between the BDS formula and the leading logarithmic approximation obtained by direct 
calculations. 

Let us now turn to the region s, S2 > 0, si , S012 1 ■S123 , S3 < where, in the QCD 
calculations, the imaginary part contains the factorization breaking term V^ut (corre- 
sponding to a BFKL ladder with the octet quantum numbers in the f-channel). In 
the BDS amplitude (174p we have, compared to the physical region with only positive 
energies, the additional phase factor C in ([75p . In the leading log approximation which 
we have described before we find (from ([74]) . ([75]) . or from (C.ll)): 

M2^4 _f\si\Y'f\s3\Y'f-iM 

r(ti)r(t3) UV UV I ^2 
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TT / / Ki2 1\ ( K23 




l-TT-^^^^-T ) • (119) 



In the imaginary part of the square brackets the e poles cancel. Comparing this result 
with (IllOp and (I105p we see that the BDS formula reproduces the lowest order term of 
the Regge cut contribution, Vcut,reduced-, but not the higher order terms (which are still 
part of the leading logarithmic approximation). We therefore conclude that, beyond 
the one loop approximation, the BDS formula does not agree with the leading log 
results listed in the previous subsection. 

Let us remark on the order 0{e) corrections in the BDS formula. As explained 
at the beginning of section 3, our analysis of the BDS formula (which applies to the 
logarithm of the scattering amplitude) does not include terms which vanish as e — > 0. 
Nevertheless, the comparison of (jllOp and (jll9p shows that such corrections cannot 
reproduce the finite (in e) terms which are missing in the BDS formula. The key point 
is that, in the BDS formula, the leading log approximation for the imaginary part of 
lnM2^4 contains terms of the order 1/e only inside 0^2 • Comparing (IllOp with ()119p 
one sees that the infrared divergent phase factor for the cut contribution in (jllOp is 
the same as in the BDS formula. Therefore, when going from lnM2_»4 to M2-^4 it is 
incorrect to expand this infrared singular piece e"*'^'^/'^, and it becomes clear that terms 
of order e in the logarithm of the scattering amplitude cannot produce constant (in e) 
terms in the scattering amplitude. As a result, our conclusion concerning the validity 
of the BDS formula is not affected by the order 0(e) corrections in the BDS formula 
for the logarithm of the scattering amplitude. 

For the 3^3 amplitude the comparison between the BDS amplitude and the high 
energy behavior in Yang Mills theories leads to the same conclusion, although some 
details are different. For the kinematic region where all energies are positive the BDS 
formula agrees with the leading log calculations, and we directly turn to the region 
s,S2 > 0, si, S3, S13, S02 < 0. The crucial element is the phase C in ([8T|) which, in 
contrast to C in ()75p . is infrared finite. Collecting, in (I80p . (I8ip . or in (D.89), within 
LLA, all imaginary parts in InMa^a we note that all terms of the form a/e cancel, and 
we arrive at: 

M3^3 ^ (\^" /^NA"' 
r(ti)r(t3) 

l + z^ f^i + a fin ^ - 1) + ..3 + a fin ^ - 1)) - z.a fin "i^'^l 1(120) 
2V \ e) \ eJJ \ (91+^3-^2)^92/^ 

The square bracket expression is infrared finite. Comparison with (jllSp shows that 
the BDS formula correctly reproduces the one loop approximation to Ucut,reduced, but 
not the higher order loops. Again, terms of order e in lnM3^3 cannot reproduce those 
finite (in e) terms which are missing in M. 

Discrepancies in the BDS finite pieces for six gluon amplitudes, starting at two 
loops, were also hinted in [39] where the equivalence between Wilson loops and MHV 
amplitudes was assumed. In a particular kinematic configuration, and for a very large 
number of external gluons at strong 't Hooft coupling, the finite pieces of the BDS 
ansatz failed when compared to the results of [50] . 
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5 Conclusions 



In this paper we have assembled the ingredients needed for the three loop corrections 
(NNLO) to the BFKL kernel in iV = 4 SYM theory at large Nc- Following earlier 
calculations we can obtain the kernel from unitarity sums, i.e. by computing squares 
of production amplitudes, keeping in mind that at large Nc the contributing diagrams 
belong to the cylinder topology. Figure [13] illustrates the production vertices which 
enter the three loop calculation. 



Regge 

s > -i 



LO 




NLO 




NNLO 





Multi- Regge 

s> IsiUsal 
k2| > -h, -t-z 







Qucisi-multi- Regge 

S> |S1|,|S3| 




Quasi-quasi- 
multi- Regge 

S > |S4| 

|si|, |S4| > |S2|, kal 

|S1|,|S4| » -t„ 




Figure 13: Diagrams contributing to the BFKL kernel in NNLLA 



Elements in the first two lines are known, whereas the building blocks in the third 
line are new: they can be (and partly have been) computed from the effective action 
summarized in section 2. 

For most of the cases, also the BDS formula can be used. The NNLO gluon trajec- 
tory function follows from the 2^2 scattering amplitude (first column); details are 
described in Appendix A. For the 2 — > 3 case (second column), the two loop approxi- 



26 



mation to the gluon production vertex can be read off from the analysis presented in 
section 3 (c.f. Eq. (1490 ). In this case one should take into account the e-corrections 
to the BDS amplitude. In column 3, we should take into account the the Reggeon + 
Reggeon 2 gluon vertex in the one loop approximation for fixed invariant masses 
of produced gluons. Based upon the analysis carried out in sections 3 and 4 we trust 
that the BDS formula for the maximal helicity violating case can be used (with e- 
corrections). For the non-maximal violating cases in column 3 one can use the results 
of [51]. Finally, in column 4 we encounter the Born vertex: Reggeon + Reggeon — > 3 
gluons for the fixed invariant mass of the gluons. This vertex has been obtained in [30] 
by means of the effective action (see also Ref. [32]). 

We have shown that the BDS amplitude M2^3 in the multi-Regge kinematics satis- 
fies the dispersive representation, which is valid in all physical regions and is compatible 
with the Steinmann relations and gluon reggeization. For the case of the gluon transi- 
tions 2 — > 4 and 3 ^ 3, in the multi-Regge kinematics and in the physical region where 
s, S2 > and si, S3 < 0, the Regge factorization of the BDS amplitude is badly violated. 
In the one loop aproximation the BDS result in this region coincides with the direct 
QCD calculations, but in higher loops we have shown that these amplitudes should 
contain the Mandelstam Regge cut in the t2-channel. It was demonstrated, that this 
cut is absent in the BDS expression and cannot be reproduced by the 0(e)-corrections 
to this expression. 

A remark is in place on the Regge-cut contribution illustrated in Fig. [12] and dis- 
cussed in section 4. In addition to the corrections to the production vertex functions 
which are illustrated in Fig. [T3] we still have to take into account those corrections 
to the production amplitude in the multi-Regge limit which do not fall into the class 
of loop corrections to the production vertices: in NLO these are just the Regge-cut 
contributions to the imaginary part in the 2 — > 4 and the 3 — > 3 cases which we have 
discussed in the previous section. The diagrams contributing to the BFKL Pomeron 
in the large Nc limit belong to the cylinder topology: two examples are illustrated in 
Fig. [m and, to begin with, we consider the discontinuity due to the 4-particle interme- 
diate state. In the left figure, on both sides of the discontinuity cut, we have the 3 — > 3 
production amplitudes continued into the physical region of a 2 — > 4 process (c.f. rhs 
of Fig. [T2j) , and in the right hand figure we recognize a configuration where the 2^4 
amplitude has to be evaluated in a region with negative energies. As discussed before, 
in the latter case the non-factorizing pieces of the 2^4 production amplitude do 
not cancel. If these contributions would survive in the total cross section, the NNLO 
BFKL Pomeron would receive a four-reggeon cut contribution, and the simple ladder 
structure would be lost. There are, however, reasons to expect that, in the large Nc 
limit, the sum of these contributions might cancel in the total cross section. Namely, 
in addition to the contribution of the 4-particle intermediate state, we also need other 
cuts, which, for example, run across one of the ladders or along one of the reggeized 
gluons above or below the cylinder. These different cuts provide similar contributions, 
but they come with different signs. It is likely that, similar to the AFS cancellation of 
Regge cuts in planar amplitudes, the four reggeon cut contributions cancel in the sum. 
We will study this in the subsequent part of our investigation. 
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Figure 14: Cylinder-type topologies in the unitarity sums for the total cross section: the 
intermediate states (discontinuity cuts) are obtained by slicing the cylinders in all possible 
ways across the intermediate momenta pa' and pb' 
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Note to be added: After our paper had been submitted, another study appeared 
which, in some parts, parahels our investigation [52]. Like ours, it studies several 
different Regge hmits of the BDS amphtudes. In the unphysical region (negative ener- 
gies), the resuhs on the energy dependence are fuhy consistent with ours. In contrast 
to our paper, however, in [52] the continuation into physical regions has not been in- 
vestigated, and the conflict with QCD calculations was not found. Recently another 
paper of these authors [57] has appeared. We completely disagree with the statement 
in their section 4.3, saying that the Steinmann relations for the BDS amplitude M5 are 
violated. In our view, the authors compute the discontinuity of the BDS amplitude 
(4.16) in an incorrect way: for example, the discontinuity in si in the physical regions 
should be defined at fixed n (this follows already from a simple one loop calculation 
in QCD where, in physical regions, the presence of simultaneous singularities in the 
overlapping channels si and S2 would contradict the gluon stability). We also find 
their section 4.6 very confusing. First, contrary to their statement, in our paper the 
dispersive representation was not used to prove the absence of the Regge factorization 
of the BDS amplitudes M2^4 and M^^^ in the physical region with s,,S2 > and 
si, S3 < 0. We simply analytically continued the BDS formula to this region and com- 
pared with the QCD results. Such a continuation is absent in the paper [5?]. Second, 
the 0{e) terms in the BDS formula do not affect any of our conclusions (this is ex- 
plained in some detail in our section 4.2). Next, the discussion around (4.29) is very 
misleading: the 'unwanted piece' in (4.29) has been derived, via the correct analytic 
continuation, from BDS, it certainly has the correct infrared properties. Comparison 
with high energy QCD calculations shows agreement with the infrared divergent piece, 
and the disagreement beyond one loop comes in when expanding, in (4.29), the finite 
term of InMg. Finally, it was shown in our Appendices C and D that the finite parts 
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of the factors C and C appear just from the analytic continuations of the dilogarithm 
functions Lz2, and they depend upon the conformal invariant cross ratios ^> and 
So what do the authors criticize? 

Recently, the paper [53] appeared where the authors calculated the three-loop Regge 
trajectory and three loop coefficient functions. Further, the breakdown of the BDS 
ansatz for the 6-point amplitude in two loops was found by direct calculations [54] in 
agreement with the predictions from the Wilson-loop calculations |55j . 

Second note to be added: Recently a new paper on the high energy behavior 
of the BDS formula appeared [58j. We do not agree with the main result of the most 
recent version 3 of this work, stating that in the multi-Regge kinematics the special 
functions appearing in the BDS ansatz are not important. To be more precise, the 
authors argue that the two limits: energy S2 ^ oo and e ^ do not commute. Ac- 
cording to appendix C, in the region s,S2 > 0, §345,5456 < the sequence of limits: 
lime_>o lims2^oo F implies that the special functions do not contribute, whereas the 
opposite order lim<j2->oo lime->o F leads to our result with the special functions being 
important. We disagree with this 'non-commutativity', since the first part of the ar- 
gument is based on a simple arithmetic mistake. Namely, starting from eq.(C.16), the 
multiplication of the factor (— P^)"*^ with (1 — $)"'^ in Eq. (C.22), in the limit S2 — > 00, 
gives the finite expression (p4_|_ +P5^)~^^, in agreement with Eq. (B.ll). In this way 
the dependence on S2 cancels out, both sequences of limits lead to the same answer 
(contrary to what is stated after Eq. (C.25)), and our result has been confirmed: in 
the BDS formula, the special functions are important in the multiregge kinematics, 
and their presence implies that the multiregge factorization is violated. 

A The 2^2 amplitude 

Let us write the BDS amplitude for the general case of n legs (see [3T]): 

00 

lnM„ = 5: a' (/«(6) (/(i)(/6) + i^«(0)) + + E^\e)) , (A.l) 
1=1 

where En^ (e) can be neglected for e ^ 0, the values of the constants are 

C7(i) = 0, (A.2) 
= -Cl/2, (A.3) 
f^'\e) = + + (A.4) 

f^' = ^g, (A.5) 
/i = -aC3/2 + a2(2C5 + 5C2C3/3), (A.6) 



7x is the cusp anomalous dimension 




and the finite remainders F^i^^ are expressed in terms of logarithms and dilogarithms. 
For the elastic scattering amplitude case we have 
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F 



(1) 



Therefore 



(1) 



2 1 

^ +ln(-t)- +ln(-s) 



In- 



+ 4C2. 



(A.c 



(A.IO) 



As a result we obtain for M4 Regge-type behaviour, as aheady discussed in the main 
part of our paper, with the gluon Regge trajectory given by 



u:{t) 



- - In ^ ) + 



1 



-t 



C2U--ln-T +• 



2e 



C3 



(A.ll) 



Note that this result at two loops is in agreement with the direct calculations [2] [3] 
based on the BFKL approach Indeed, in Ref. [3j the following expression for the 
gluon Regge trajectory was obtained in the MS-scheme (using the same notations): 



In ^ I + a 



^-C2)(--ln^)+^-| 



6 



2e 



9 



(A.12) 



The contribution of the scalar loop to this trajectory is proportional to the contribution 
of the fermion loop [3] 



LO- 



US 



it) 



1 



1 



^-ln^-- - + 21n— T 
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9 



(A.13) 



where Us is the number of scalar fields transforming according to the adjoint represen- 
tation of the gauge group. For the transition from the MS-scheme to the dimensional 
reduction (DRED) scheme, which respects = 4 supersymmetry, one should first 
increase the number of scalar fields 



6 + 2e, 



(A.14) 



because, in the pure gluonic contribution. An = —2e for the gluon fields was taken 
into account after performing the dimensional regularization 4 — > 4 — 2e. This gives 
the additional contribution to ujy^{t) 



Au 



MS 



(A.15) 



After that the subsequent finite renormalization of the coupling constant needed for 
the transition between the MS and DRED schemes 



1 2 
a a 

6 



leads to the above result for the trajectory 



(t) ^ u;it) 



In- 



<2 



1 

2^ 



In- 



2 



(A.16) 



(A.17) 



Concerning the residues T{t) of the Regge pole, they have been calculated in the 
one-loop approximation in QCD [33]. In supersymmetric models the helicity non- 
conserving contribution of each of the colliding gluons is cancelled, in accordance with 
the BDS ansatz. 
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B The 2^3 amplitude 

For the 2 — > 3 production amplitude we have (see Figl3|) 

f(l). ^ 5 1 , i-s){-Sl)i-S2)i-h){-t2) 



2e2 2e /iio 

- lfln2Z£+in2Z£l+in2Z£^+ln2Z| + in2Z|l (B.l) 
4 V /^^ Ai^ fJ-^ /"^ fJ.^ J 

i^i^^ = — -In In — -In — - In — -In In ■ ^ 
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4 — si — S2 4 —ti —si 4 — s — ti 
iln^ln^-iln^ln^ + i^C2. (B.2) 



Thus the total contribution in multi-Regge kinematics is 

5 



2e2 



_ i Zii + 1 ,„ fi,, (zMzM _ i) 1 Zjl + 15c, . (B,3) 

4 /i2 2 /i2 V e y 4 4 

In this way we obtain the Regge factorization of the production amplitudes, discussed in 
the main text. Let us note that, formally, this result is exact and the amplitude can be 
written in this factorized form in all five channels obtained by the cyclic transmutation 
of the invariants s, ti, si, S2, ^2- 

In the one-loop approximation in QCD the Reggeon-Reggeon-gluon vertex contains, 
apart from the Born structure proportional to the vector C{q2,qi), also the contribution 
proportional to the gauge-invariant vector ^ — ^ [33j . In the supersymmetric theories 
this contribution is cancelled, in agreement with the BDS ansatz. 

One can calculate also the production amplitude in the quasi-elastic kinematics, 
where s ~ si ^ S2 ~ ti,t2, k'^- The amplitude here has the usual Regge factorization. 



C The 2^4 amphtude 

In the case of the 2 — > 4 transition we have (see Fig. [6]) 

f(l). ^ _ 3 1 {-s){-Sl){-S2){-S3){-h){-t^) 

1 -S , ,2 , i„2 -S2 , ,2 , ,2 , i„2 -^3 



^ , In^ ^ + In" ^ + In" ^ + In" ^ + In" ^ + In" (C.l) 

4 V ^J'^ /^^ /"^ ' 

Fr^^ = In In In — - In — In — In — 



2 — S012 — S012 2 —t2 —t2 2 — Si23 —^123 
1 1 -52 , -Si 1 . -Si , -h 1 -tl -S 

— m in In in in m ■ 



2 — S012 — S012 2 —t2 —t2 2 — Si23 — Sl23 

- ^-U2(l-^^\--U2(l-^\--U2(l-^ 
2 V SOI2S123/ 2 V *2S012/ 2 V i2Sl23 

4V -8012/ 4V -5123/ 4V -5012/ 2^ ' 
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where the dilogarithm function is defined as 

df 

Li2{z) = - -ln(l-t). (C.3) 
Jo t 

In multi-Regge kinematics it is natural to introduce the independent variables 

Si, S2, S3 , -Ki2 = —7 ^ , -K23 = ^ , ^ = 7 77 7 • (C.4) 

(-SOI2) (-S123) (-S012)(-Sl23) 

Note that the variable $ is unity in the region where all above invariants arc negative, 
but ^> = cxp(— 27r7;) in the physical region where s, S2 > 0, Soi2jSi23 < 0. In the 
multi-Regge kinematics we obtain the following general result: 

+ In — ^ (- - In + In (- - In — + In — ^ (- - In 



-[l^^Z^+ In^ - In ^ fin iz^l^zM - 1 U 1 In 



2 /x2 V e; 4 V -t2 -t2) 2^' 



At first sight the arguments of the dilogarithm functions in the multi-Regge kine- 
matics are either or 1, and we can use the relations 

Lz2(0) = 0, Li2(l) = C2. (C.6) 

However, in the physical region s,S2 > 0, si, sa, soi2) si23 < it is needed to be 
cautious: we should analytically continue the expression 

/($) = Lz2(l = (C.7) 

S012S123 

from the region $ 1 to the region $ ^ e~^'^' along a unit circle. In multi-Regge 
kinematics we have 

_ Soi2'Sl23 /r , \2 /r^ o\ 

S2 ~ [kl + /C2 j (C.8) 

and 

S2(i = +^2) ^ (C.9) 

Therefore 1 — $ > 0, and after the analytic continuation we obtain 

/($) = -/ ^ln(l-t) + 2Tri -f«27riln(l-^>) (C.IO) 

Jo t Jit 

with ln(l — $) being real valued. We obtain the following result in the physical region 

S,S2 > 0, Si,S3,Soi2,Sl23 < 0: 



{-ti){-U) 1 



-'6 (ej + i^e - -^ + ^^ 

+ ln^fi-ln4Uln 

+ - In — ^ ( In — + In — ^ ) + - In — ^ ( In — ^ + In — ^ 

,2-tl , ,„2-t3\ , 1 ,„ (-tl)(-t3) , 7, 
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It is possible to derive, from the BDS amplitude, an expression for M2_»4 in the 
one-loop approximation for the quasi-multi-Regge kinematics, where s >> si,S3 >> 
S2 ~ ^1,^25^3- In this case it is convenient to introduce Sudakov variables for the 
momenta of the two produced particles 

kr = PrPA + OirPB + , (kr f = -kl, (C.12) 

where 

„2 2 

1 > /?i ~ /?2 > — , 1 > ai ~ 02 » — , sarl3r = kl ^ qlr^ ql^ qlr^ . (C.13) 
s s 

We can express various invariants in terms of these variables 

S2 ~ s{(3i + (32){ai + 132) - ih + k2? , (C.14) 

si « ais , S3 f« (32S , S012 ~ (ai + "2)5 , S123 ~ (^1 + /32)s • (C.15) 

The expression for the function /(e) = /6^''(e) + Fq^^ (for <I> = 1) in the quasi-multi- 
Regge kinematics can be obtained by adding an aditional term 



f{e)^ f{e) + Af, (C.16) 



where 



A/ - 1 SK12K23 stit^, ^ 1 S012K12 hsiS2S 



2 S1S2S3 S012Sl23/^ 2 S1S2 i2Soi2'5l23 

, 1 , Sl23/«23 , ilS3S2S 1 , 2 S«^12K23 1 , 2 •S012K12 1 , 2 ■5l23«^23 

H — In in 5 in m In 

2 S3S2 t2S012Si23 4 S1S2S3 2 S1S2 2 S3S2 



1 SK12K23 S012'Sl23^12'^23 1 , S012'^12 , Sl23'^23 

- — m m 7) m m 

2 S1S2S3 SiS3S^ 2 S1S2 S3'52 



-h^2 ( 1 - ) - h^2 ( 1 - - h^2 ( 1 - . (C.17) 

2 V SOI2S123/ 2 V ^25012/ 2 V t2Sl23/ 

Here the signs —1 are implied to be in front of all invariants Sj,tj. Note that the ex- 
pression for A/ in the quasi-multi-Regge kinematics does not contain large logarithms, 
because the arguments of all logarithms and dilogarithm functions are of the order of 
unity. It is proportional to the logarithm of the amplitude for the transition of two 
Reggeized gluons into two particles with the same helicity. Similar to the case of M4 
and M5 the expression for M2^4 in the quasi-multi-Regge kinematics coincides with 
the exact BDS amplitude. The transition of two reggeons to particles with opposite 
helicity in the one-loop approximation can be found in Ref. [51]. These transition 
amplitudes are needed for the calculation of the next-to-next-to leading corrections to 
the BFKL equation. 

D The 3^3 amplitude 

Here we consider the BDS amplitude Mq in the channel corresponding to the transition 
3 — > 3 with the following invariants (see Fig. [9l): 

s = ipA + ki+ psf , si = {pa + kif , S3 = (pb' + k2f , (D.l) 

Sl3 = (^1 + Psf , Sq2 = ipA' + k2f , t2 = ipA' + k2 - PA^ , (D-2) 
h = {PA' - PAf , ^3 = {PB' - Psf , t2 = {pA' -PA- kif . (D.3) 
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The functions and Fq^^ in this case are given by |31j : 

f{l)fs _ 3 1 i„ i-Sl)i-Sl3){-S3){-S02){-ti){-t3) 



_l/'l„^Z^l + l„2Z£H + l„2Z|L + l„2Z^ + ,„2Z^i + ta2Z|| , (D.4) 



/i^ ^2 ^2 ^2 ^2 ^2 

Fg = — — In i In — — — In — ^ In — 4- — — In — — In ■ ^ 



-lln 

2 


— s 


-Sl3 
— S 


-il„ 

2 




In - 


S02 





^2 ~^2 ^ "^2 
•502 1 1 , 1 -•51 

in — ; In in 



s -s 2 -t'2 -t'2 2 -t2 -t2 

In multi-Regge kinematics 

- s > -si, -S3, -t'2 > -ti, -t2, -is > (D.6) 
it is helpful to use the definitions 

{-Sl){-t'2) ^ {-S3){-t'2) w (-^13)(-^02) 

-•502 -^513 (-*2)(-s) 

which allows us to simplify the above expressions 



-1 fin^ ^ + - In ^ (m th}(-hl - iU - m ^'''^^-''^ 



2 V /"^ ey 4 V -t2 -^27 2^' 



In the physical region, where s,i2 > and si, S3, S02, •513 < 0, one has = exp(27ri), 
i.e. we have to continue in along the unit circle. The relation 

- ^'2(1 - « (gi + 93 - 92)2 p_g) 

implies that, after continuation, 1 — < 0. Therefore, Li2 becomes 

A^^*' (it /-i-*' ^+ 

/($') = -/ _ in(l - t) - 27ri / — « 27r2 - 27riln |1 - $'| , (D.IO) 
Jo t Ji t 

which allows to obtain the extra phase factor C' violating the Regge factorization in 
this physical region. 



E High energy scattering amplitudes in the 
leading logarithmic approximation 

In this appendix we briefly summarize results for the high energy 2 ^ 3, 2 — > 4, 
and 3 — > 3 scattering amplitudes in Yang-Mills theories in the leading logarithmic 
approximation. 
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For the 2^3 case most of the results have already been listed in section 4.2. We 
only quote, for the physical region where all energies are positive, the factorized form: 

^2^3 2S 



sin7r(a'i — UJ2) 

Here the trajectory functions uJi and the production vertices Vi, V2 have been com- 
puted in LLA and NLO, whereas the phases and sin factors are part of the analytic 
representation, and do not need to be expanded in powers of g^. However, since in this 
paper we restrict ourselves to the LLA, we can put (e~*'^Ki2)~'^^ ~ 1. For the real part 
we find the result (j96p which coincides with (4). 

For the 2^4 amplitude we start from the ansatz ([9] 

^2^4 2s f dio'2 



r(ti)r(t3) tit2t3J 27ri 

5^|\'^l-'-'2 /|sp^2|\'^2-'^3 "^3 Wi{ti,t2,t3,Ki2,K23;i^2) 

';i2'?2'3?3 



/i^ / V /^^ / V/U^y " " sm 7rwi2' sm 7ra;2'3 

Ny"^ /Nldy^""' f ^ V' 62'e2'iei ^2(tl,t2,^3,^12,^23;c^2) 

/i^ / \ fJ''^ J \fJ''^J sin7rti;32' sin7rix'2'i 

821^''"^' fMY'-'^' f\s\Y\ VF3(tl,t2,t3,«12,«23,'^123;c^^) 

42'ia3?3 



/i^ / V /U^ / KfJ''^/ ' ' ' sm7ru;2'i sm7ra;i3 

V /i^ / \ fJ-"^ J \^Jp') sin7ru;2'3 sin7rw3i 

V /i^ / V A* / VAi / sm 7ra;32' sm 7ru;i2' J 

The partial wave functions VFi=i,2,3,4,5 have been listed in section 4.2. They have 
been obtained from the five single energy discontinuity equations, and use has been 
made of the BFKL bootstrap equations in the color octet channel. 

Inserting these partial waves into the ansatz (|97p or (jE.2p . we can study the full 
amplitude in the different kinematic regions. From now on we will specialize on the 
planar approximation, i.e. in the signature factors in eq. (|93p we only retain the phases. 

Beginning with the physical region where all energies are positive, we first collect 
the Regge pole terms in all five partial waves Wj. Their sum can be written in the 
simple factorizing form: 

^2^A,pole _ 2s 
r(ti)r(t3) ~ tit2t3' 

^C-^^\,^\Y2 (e-'"^12)-"^^l(^l,t2,^12)-(e--"^12)'"^1^2(tl,t2,ACi2) , „ _ 

sin 7ra;i2 

(e'-"/^23)-"^V^lfe,t3,At32)-(e---^23)-"^^2(t2,t3,K32) ^^„,^j„^j^^3^^^^3^ 

sin 7ra;23 

In order to arrive at this result, we have combined, in (jlOOp and (jlOip . the Regge pole 
contributions of and W4 (together with the signature factors), and we have used 
the identity: 

sin7ra;23 sinvraii ^ sin7r6<j2i sin7ru;3 ^ (E 4) 

sin7ra;i3 sin7ra;2 sin7rti;3i sin7ra;2 
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The production vertices are the same as in the 2 — > 3 case, (jE.ip . As in the 2 — > 3 case, 
in the leading order approximation we put, in ()E.3p . {e~™ k)~'^ ~ 1. For the real part 
the terms proportional to ln(«;/;U^) — ^ in the production vertices Vi and V2 cancel, 
and we are back to the factorizing form in (4). As to the addititional Regge cut pieces 
contained in and they cancel completely: 



A2- 



*4,cut 



r(ii)r(t3) 

-2s f duj'n , ,x,„' / 1 



(Niir / |^(e-^is2ir^ + {Vcut - V,) (|s3ir (e.5) 

J Zm VsmTroiis smTrwsi/ 



0. 



It is instructive to study the cancellation of the imaginary part of this Regge cut piece 
in more detail: from the representation (j97p which shows the energy phase factors 
explicitly it is straightforward to compute the single discontinuities in S2, soi2) S123, 
and in s. When summing these single discontinuities (i.e. when computing the full 
imaginary part), we find complete cancellation of the Regge cut piece. This cancellation 
of Regge cut contributions in the planar amplitude is nothing else but the Mandelstam 
mechanism |45j of the cancellation of the Amati-Fubini-Stanghellini Regge cut [M] 
in planar diagrams. 

The unphysical region where all energies are negative can be obtained from ()E.3|) 
and (jE.SP by simply putting the phase factors equal to unity: the factorizing form of 
the Regge pole contributions is preserved, and the cut pieces in W3 and VF4 cancel. 

Most interesting is the physical region where s,S2 > and si, S3, soi2) ■S123 < 0. 
Nonzero phases appear only in s and in S2- After some algebra we rewrite ()E.2[1 in the 
following form: 

-(isiircNsir- (E.6) 



e-'^'^H\s2\) 



U12 



r(ti)r(t3) tit2t3 



_^g-.«a.i2_ z g-./.a.32_ i ^ 



sin7rti;i2 sm7ru;2i/ \ sm7ru;23 sm7ru;32 



\ V2V1 



sm 1^02 sm 7ra;i3 sm 11002 sm ttuj^i J sm 110021 sm 11002^ 

2s „ f dojL /e-^'^^Sl e-i7ra;i3 ' 



^(e— |s2|r^ _ + _ _ v^)) . 

iTTi y sm7ra;i3 sm7ra;3i I 

It is important to study the infrared singularities of the phases of this expression. First 
we note that the prefactor e"™'^'^ contains, in lowest order in a, the 1/e singularity of 
the gluon trajectory (j44p . All other phase factors contain differences of trajectory 
functions and are finite as e — > 0. Expanding the square brackets in powers of a we 
arrive at (|107p : 

Vcfe,gi)c7fe,g2)(Niir(N3ir- 



r(ti)r(t3) tit2h^ 



e-*""2(|s2|r 



vr / K12 
1 + \ uJi +U2 + a ^In — 2" 



1\ A «;23 1' 
- + (J3 + W2 + a In 



-2iTT- 



2s 



^1*2*3 
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Following the steps described in section 4.2 we factor out the gluon trajectory (details 
are presented in [i7]): 



^{e~^-\s2\r2V,^t = 9^C{q2,qi)C{qs,q2){e-'^\s2\r J ^{e-'^\s2\r''V,ut, reduced- 

(E.8) 

Insertion into ()E.7p leads: 



.vr 



1 + i— { LUi + UJ2 + a { In 



Kl2 



1 

e 

-2i7r 



+ 0^3 + c<j2 + a ( In ^ 



dui2 
27ri 



1 
e 

reduced 



(E.c 



The one loop approximation of the Regge cut contribution, (jlOSp . is infrared singular: 



9^C{q2,qi)C{q3,q2) / ^(e '''|s2ir2T4„^^^^^„^^^ 
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2C{q2,qi)Ciq3,q2) 



a ( In ^^^^f^ + - 
V mi e 



+ ... 



(E.IO) 



whereas the two loop and higher order terms of Vcut,reduced can been shown to be 
infrared finite [47]. As a consequence, in (jE.Oh the coefficient in the square brackets is 
infrared finite, and the singularities are collected in the overall phase factor e"*'^^^. 

For completeness we also list those energy discontinuities which do not vanish in 
this kinematic region, the discontinuity in the total energy s and the discontinuity in 
S2- We again start from the analytic representation (|E.2[1 . After some algebra (which 
includes approximating phase factors by unity) we finqj: 



disc. 



A 



2^A,cut 



-ITS 



r(ti)r(t3) tit2t3 



i\si\T 



du2 
27ri 



i\s2\r^ Vcut {\s3\r. 



In a similar way we compute the discontinuity in S2'- 



discc, 



A 



2^i,cut 



-ITS 



T(ti)r(t3) ht2t3 

where, instead of (|103p . 

d^kd^k' 



^(k2i)-^ v^ut {\s3\r, 



(E.ll) 



(E.12) 



Vr 



cut 



t2Nc 4 



[C{q2,qi] 



ql 



(27r)6 ^-^-'-^ ^k + qrf 
X G'^^^\k + q2,-k;k' +q2,-k'-J2) 



C{k + q2,k + qi) 



X [C{q3,q2)-C{k' + q3,k' + q2) 



4 



{k' + q3f 



(E.13) 



In lowest order we have: 



T>(0) 2C{q2,qi)C{q3,q2) (, Ki23/i^ , 1 

^cut = 9 — a I In + 



l^l2l<-23 



*We use the definition discf{s) — jr (/(s + ie) — f{s ~ ie)). 



(E.14) 
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A completely analogous discussion applies to the case 3 — > 3 (Figs. [U [TT]) in the 
multi-Regge region (for a detailed discussion of the 'analytic' representation see IMj)- 
Our ansatz is: 

^3->3 2s f dw'2 



/i^ y \ iJ? ) V/^^/ ' ' sm 7rwi2' sm 7rLiJ2'3 

|S3| V"V|S2|\"^""VI^I^"\ t C f/2(tl,t2,t3,'^12,At23;^) 

-T ";;r 7:2 ) ?32'42'i4i 



IJL^ J \ 1-^ J XtJ- / sin7rw32' sin7ra;2'i 

^ ^ * - J2) 



I |\ LuL—UJ-A /I |\ iuL—UJl /I |\ LUl-^UJ-A—iuL T T f J. J. J. /~ 

|S02| \ ^ /I'^ish Vl'^h V ^ ^ (^3(il,i2,t3, ^^12, ^23, '«123;'^2. 



2 / \ 2 ] \ 2 I 6'l6'3C(l+3)2' . . . / , 

/ V / V/^ / sin7rw2'i sin7rw2'3 sin7r(wi + 0^3 — 0^2' 

|S02|\"^*^^ /kl3|\"^*^^ ^k2|\"^-"^-"V , , f/4(tl,i2,*3,Atl2,/^23,Kl23;^^) 



66'^2'(l+3) 



/ \ IJ? J V A*^ / - - - K — . Sin7r(^LiJ2' — Wi — W3j 

|g3|\ "^s-i^a '^i-'-'2 '-'2 C/5(il,t2,t3,Kl2,K23;'^2)l lT?^r,\ 

-9 ?32'?12'?2' : : , (ll..i5j 

where 

= e— + 1, = e-''^«'^»+'^^)-'^'=) + 1. (E.16) 

The pieces labelled by 1,2,5 are 'normal' and contain only Regge poles. They coincide 
with those of the 2 — > 4 amplitude: 

Ui = Wi, i = 1,2,5, (E.17) 

and they fit into the factorization pattern, The terms 3 and 4 have the extra Regge 
cut piece shown in Fig. [12] (right figure), which is described in terms of the color octet 
BFKL equation. In analogy with (jlOOp . (|10ip one finds: 

sin vrwi sin vrtJs _^ ^ , , 1 , , 

U3 = : V2{ti,t2,Ki2)—, Vi{t2,t3,K23) 

Sm TTUJ2 ^2 ~ "^2 

+ sin7ra;2'i sin7ru;2'3 (f^cut - Up) , (E.18) 
1 1 

Ua = T^2{tl,t2,K.l2)— Vi(t2,t3, «^23) 

sm t:oj2 UJ2 ~ "^2 

+ (Ucnt-Up). (E.19) 

The cut piece has been given in pU3|) . (|1U5|) ): 

t2N, 4 f (fkd?k' ql 
= J (2^)6 (fc-,iP^^^^-^''^^-^) 

X G^^^\k,q2-k;k',q2-k';J2)C{k' -k2,k')j^^^, (E.20) 
with the lowest order approximation 

rr(0) _ g^g(g2,gl)C(g3,g2) , ^12^23 ,^ ^ 

f^cut - IT-, am——— ^.2 2 - (l^-21j 

2^2 [qi +qz- q2rq2 

It contains Regge cut singularitities, and breaks the factorization. Note that, in con- 
trast to the 2 — > 4 case, the one loop approximation of the Regge cut term, has no 1/e 
pole, i.e. it is infrared finite. 
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In analogy with the 2 — > 4 case, these Regge cut pieces does not show up in the 
physical region where all energies are positive. It is, again, only in the other physical 
region s,S2 > 0, si, S3, S13, so2 < where these pieces become visible. Proceeding in 
the same fashion as before (1E.6P we find for this region: 



(i^2ir 



{\si\ri\s3\r 



Sm TTLU12 



+ e 



2i 



V1V2 



sm 7ra;2 

The last term can be written as 



+ 



sm 7raJ2i 
2is 



sm 7rci;23 sm vrci;32 



^1*2^3 



(i«iir(N3ir 



2TTi 



ie~'^\s2\r^U,^t = g'Ciq2, qi)Ciq3, g2)(e-*"|s2| 



duj2 
27ri 

dw'2 
27d 



It is important to note that the one loop approximation, U^^^, 

dujn 



(e-^-|s2|r2[7,„t.(E.22) 



(e ^'^\s2\Y'^Ucut,reduced- 

(E.23) 
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g'C{q2,qi)C{qs,q2) f ^(e-'^ksir^f/, 
,2C{q2,qi)C{q3,q2) 



o In ■ 



'^12^23 



cut,reduced 



(E.24) 



{qi +93 - 92)^^2' 

as well as the higher order terms are infrared finite. We therefore write ^3^3 in the 
following form: 

^3^3 2s 



r(ii)r(t3) ^1*2*3 



i\si\n\s2\r'i\s3\r9'C{q2,qi)Ciqs,q2 



.TT 



KI2 



1 



1 + z— LJi + a In + W3 + a In 



K23 



1 



duj2 
21x1 



\^ \^2\) '~' cut, reduced 



(E.25) 



On the rhs, the square bracket term is infrared finite. This shows that the infrared 
structure of ^3^3 is quite different from ^2^4. 

We conclude this section by listing the discontinuities in the energies s and S2: 



disc. 



A' 



3^3, cut 



-vrs 



r(ti)r(t3) tit2h 



{\si\r 



dw2 
27ri 



{\s2\r- Ucut {\S3\) 



1^3 



In a similar way we compute the discontinuity in S2: 



disCs 



A 



3^3, cut 



-vrs 



r(ti)r(t3) tit2t3 



i\si\) 



di / 

^(k2ir^ Ucut {\s3\r, 



(E.26) 



(E.27) 



where U is obtained from U in the same way as V was obtained from V. 

In order to compare these results with the BDS formula, we divide the scattering 
amplitudes by their Born approximation. For example, we obtain M2_+3 by dividing 
^2^3 by the Born approximation 

2s 

—9^SxA\A'^^B\B'9C{q2,qi) 

tlt2 

and M2^4 by dividing ^2-*4 by 



2s 



^1*2^3 



5^<5AAA4,'^AsAB,5^C'(g2,gi)C(g3,g2)- 
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